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Q ■ We revisit the earlier determination of as{Mz) via perturbative analyses of short- 

distance-sensitive lattice observables, incorporating new lattice data and perform- 
ing a modified version of the original analysis. We focus on two high-intrinsic- 
scale observables, log(Tyii) and log(Tyi2), and one lower-intrinsic-scale observable, 
j^ ' log(M^i2/no), finding improved consistency among the values extracted using the dif- 

ferent observables and a final result, as{Mz) = 0.1192 it 0.0011, ~ 2a higher than 
the earlier result, in excellent agreement with recent non-lattice determinations and, 

,^ ] in addition, in good agreement with the results of a similar, but not identical, re- 

analysis by the HPQCD collaboration. A discussion of the relation between the two 

cn ■ re-analyses is given, focussing on the complementary aspects of the two approaches. 
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PACS numbers: 11.15.Ha,12.38.Aw,12.38.Gc 

I. INTRODUCTION 

The strong coupling a^ is usually characterized by giving the value, as{Mz), in the 



k>( I MS scheme at the conventionally chosen Uf = 5 reference scale /x = Mz- A high precision 

^ ■ determination of as{Mz) based on the perturbative analysis of short-distance-sensitive 

lattice observables computed using the a ~ 0.09, 0.12 and 0.18 fm n/ = 2 + 1 MILC data 
was presented in Ref. [1]. The result, as{Mz) = 0.1170(12), plays a dominant role in 
fixing the central value of the current PDG assessment [2], as{Mz) = 0.1176(20). 

Over the last year, a number of improved non-lattice determinations of C(J.^zl have 
app eared, in a variety of independent processes, over a wide range of scales [4, |5|, la, [Zl, M, M, 
10|, lul, 112, ll3l] . The results, given in Tabled (with all errors combined in quadrature), yield 



a weighted average, as{Mz) = 0.1190(10), ~ 2a higher than the lattice determination. 
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TABLE I: Recent non-lattice determinations of as{Mz] 



Source 



asiMz) 



Global EW fit [4, 5} 

Hl+ZEUS NLO inclusive jets [6] 

HI higli-Q2 NLO jets [7] 

NNLO LEP event shapes [8] 

NNNLL ALEPH+OPAL thrust distributions [9j 0.1172 ± 0.0022 

a[e+e- -^ hadrons] (2-10.6 GeV) [10] 



0.1191 ±0.0027 
0.1198 ±0.0032 
0.1182 ±0.0045 
0.1240 ±0.0033 



r[T(i£)-nX] 
r[T(is)^x] 



m 



hadronic r decay [12, 12, 1^ 



iiQn+O-0090 
iiqn+0.0060 

0.1187 ±0.0016 



This difference, though not large, motivates revisiting the lattice analysis, especially in 
light of the existence of new high-scale (a ~ 0.06 fm) lattice data not available at the 
time of the earlier study. We perform such an extended re-analysis in this paper. 

The rest of the paper is organized as follows. In Section [TTl we outline the original 
analysis, specify our own strategy for implementing the underlying approach, and clar- 
ify the difference between our implementation and that of the earlier study and recent 
HPQCD re-analysis. In Section IIIIl we discuss the details of, and input to, our version 
of the analysis. Finally, in Section HVl we present and discuss our results. 

II. THE LATTICE DETERMINATION OF as[Mz) 
A. The original HPQCD/UKQCD analysis 



In Ref. [ij, as{Mz) was extracted by studying perturbative expansions for a number 
of UV-sensitive lattice observables. Ok- The generic form of this expansion is 



Ofc = ^ ciJ^ay(Qfe)^ = DkaviQk 






N=l 



AI=0 



aviQkf 



where Qk = dk/a are the Brodsky-Lepage-Mackenzie (BLM) scales p^] for the 0^, and 



.(fc) 



= 1. The coefficients c]^2 3 (equrvalently, Dk,Ci , 



and C2 ) have been computed in 



3-loop lattice perturbation theor y |16l| . and, with the corresponding dk, tabulated for 
a number of Ok in Refs. [l|, [ly, llTI]. In Eq. ([1]), ay(/i) is a coupling with the same 
expansion to O(a^) (with as the MS coupling) as the heavy quark potential coupling, 
, but differing from it, beginning at O(a^), in a way that will be specified below. The 

and hence the (3 function of ay, defined in 
^Q/3^ay'^^(/i), with av = ay/vr, is determined 



a 



y 



expansion coefficients are known to 0{al) 

our conventions by yU^(iav(/i)/(i/U^ = —J2n 

to 4 loops by the known coefficients, /?o 

coefficients c[ , C2 , and C3 

Ok in terms of any variable, a^, sharing the same expansion as ay out to 0{a 



/?3, of the 4-loop MS (3 function [18|. The 



and &^' tabulated in Refs. [l|, |l6|, ll7| are valid for expansions of the 



With only the known, third order terms in the expansions of the Ok-, no value for the 
reference scale coupling, ayl^-S GeV) = a^, was found to produce a simultaneous fit to 
the data at all three lattice spacings employed [ij. In consequence, terms out to tenth 
order in the expansion of Eq. ([1]) were incorporated, the unknown coefficients ci ... ^g being 
fitted using input Bayesian prior constraints. The 4-loop version of /3^ was used to run 
ay to the scales Qk relevant to each of the given observables at each of the three lattice 
spacings. Linear extrapolation in the quark masses was employed, and possible residual 
mass-independent non-perturbative (NP) contributions estimated, and subtracted, using 
the known leading-order gluon condensate contributions to the relevant Wilson loops [19i] . 

The scales ri/a and ri, which determine the lattice spacing, a, in physical units, as 
well as the gluon condensate, {asG^/n), required for the mass-independent NP subtrac- 
tion, were determined as part of the independent fit performed for each of the Ok- This 
was accomplished using an augmented x^ function in which the squared deviations of 
the relevant parameters from their input central values were scaled by the squares of the 
input prior widths. For ri/a and ri the central values and widths were provided by the 
measured values and their uncertainties. For (a^G^/Tr), a central value and uncertainty 
±0.010 GeV^ (~ the conventional SVZ value 0.012 GeV^ 0) were employed |2l|. While 
this procedure allows ri/a and ri (which should be characteristic of the lattice under con- 
sideration) to take on values which vary slightly with the Ok being analyzed, one should 
bear in mind that the measured uncertainties, which set the range of these variations, are 
small compared to the variation of scales across the a ~ 0.09,0.12 and 0.18 fm lattices 
employed in the analysis. The impact of any potential unphysical observable-dependence 
of the physical scales on the fitted ay and Cn should thus be safely negligible. The 
situation with regard to the independent fitting of {asC^/ir) for each Ok is potentially 
more complicated, and will be discussed further below. 

The resulting best fit value for ay, averaged over the various observables, was then 
matched to the n/ = 3 MS coupling, and the corresponding Uf = 5 result, as{Mz), 
obtained via standard running and matching at the fiavor thresholds [22, l23|, yielding 
the result, as{Mz) = 0.1170(12), already quoted above. 

Regarding the conversion from ay to ag, one should bear in mind that, while the 
expansion for ay in terms of a^ is, in principle, defined to all orders (see below for 
more on this point), the coefficients beyond 0{af) involve the currently unknown MS (3 
function coefficients /34, /^s, ■ ■ ■. The nj = 3 conversion step is thus subject to a (hopefully 
small) higher order perturbative uncertainty. As will be explained in Section III CI with 
the definition of ay employed in Ref. [l!], the higher order perturbative uncertainties are, 
in fact, entirely isolated in the V — > MS conversion step of the analysis. 



B. An alternate implementation of the HPQCD/UKQCD approach 

The higher order perturbative uncertainty encountered in matching ay to as can be 
removed entirely by working with any expansion parameter, ax-, whose expansion in a^ 
is fully specified. We take ar to be defined by the third-order-truncated form of the 



relation between ay{fi^) and as(yU^) [2^ which, for Uf = 3, yields 

arifx^) = a,(/i2) [l + 0. 5570a, (fx^) + 1.702a^^{fj.^)] . (2) 

The jS function for a^, (3"^, is then determined to 4-loops by the known values of 
Po, ■ ■ ■ , P3. With all coefficients on the RHS positive, ax runs much faster than a,, a 
fact reflected in the significantly larger values of the non-universal f3 function coefficients, 
j3j = 33.969 and Pj = —324.393. This makes running ar using the 4-loop-truncated /?^ 
function typically unreliable at the BLM scales corresponding to the coarsest (a ~ 0.18 
fm) lattices considered here. Since, however, the 4-loop-truncated MS running of a, 
remains reliable down to these scales, and the relation, Eq. ([2]) is, by definition, exact, 
the running of ot may be performed by converting from ot to a, at the initial scale, 
running a, to the final scale, and then converting back to ax- This proceedure will be 
especially reliable for )k like log(iyii) and log(Vri2) with lowest BLM scales > 3 GeV. 

Though the conversion from the fitted reference scale ax value to the equivalent MS 
coupling as can be accomplished without perturbative uncertainties, higher order pertur- 
bative uncertainties do remain in the analysis. To see where, define ao = ctriQa), with 
Qo = Qfc*°^ = dk/ttmin the maximum of the BLM scales (corresponding to the finest of the 
lattice spacings, amin) for the observable in question. Expanding the couplings at those 
BLM scales corresponding to coarser lattices, but the same observable, in the standard 
manner as a power series in ao, ariQk) = '^N=iP'^(^k)c(o (where tk = log (Ql/Ql), and 
the PAr(t) are polynomials in t), one finds, on substitution into Eq. ([T]), 



Dk 



l^k 



■ + a^ (cf + ■■■)+ a^ (cf - 2.87cPtk + ■■■)+ a^ (cf - 0.0033/3^^, 

-3.58cfhk + [5.13tl - 1.62tfc]cf ^ + ■■■)+ al U^^'^ - 0.0010/?Jtfe 

+ [0.0094t^ - 0m65cfhk]Pj - 4.30cf ^tfc + [7.69t^ - 2.034]cf ^ 

+ [-7.354 + 6.39t^ - 4.38tfc]cf ^ + •■■)+•■■. (3) 

where the known numerical values of Pq, ■ ■ ■ , /3j have been employed, and we display 

only terms involving one or more of the unknown quantities /?J, /^l", ■ ■ ■, Cg , C4 , • • •. 

Running the MS coupling numerically using the 4-loop-truncated (3 function is equiv- 
alent to keeping terms involving Pq,- ■ ■ , jS^ to all orders, and setting Pi = P5 = ■ ■ ■ = 0. 
The neglect of P^, P^- ■ ■ means that /3f , /3j, • • • also do not take on their correct physical 
values, leading to an alteration of the true t^-dependence, beginning at ©(ag). Since it 
is the scale-dependence of Ok which is used to fit the unknown coefficients C34 ..., as well 
as ao, we see immediately that the 4- loop truncation necessarily forces compensating 
changes in at least the coefficients C45.... A shift in q , however, also alters the 0(ao) 

coefficient, which will, in general, necessitate an approximate compensating shift in c\ 
as well, and, in consequence, a further compensating shift in ao. From Eq. ([3]), the size 
of such effects, associated with the truncation of the running, and unavoidable at some 
level, can be minimized by taking Qo as large as possible (achieved by working with the 
observable with the highest intrinsic BLM scale) and keeping tk from becoming too large 
(achieved by restricting one's attention, if possible, to a subset of finer lattices) [25i] . 



C. More on the relation between the two implementations 



For Uf = 3, in our notation, the relation between ay and as, to O(a^), is [2J 

^■^q') = asifi') [I + ^,{fiyq')as{fi') + ^2{fiVq')as{fi')] (4) 



«yv 



where Kalx) = [as + 1Q(3^ log^(a;) + (16/?i + Spoai) log(x)] /levr^, with 03 = ^ + Sevr^ - 
Ivr"^ + 14^(3), and ki{x) = [7 + 4/3olog(x)] /47r. Our expansion parameter, 0^(9^) is 
defined to be equal to the RHS of Eq. (jl]) with jj? = q^, leading to the numerical result 
given in Eq. ([2]). The conversion from ut to as can be performed exactly but the absence 
in j3j^ ... of terms ex /34,5,... induces a perturbative uncertainty in the values of our fitted 
parameters, one which can, however, be reduced by working with high scale observables 
and fine lattices. It is also possible to test for its presence by expanding the fits to include 
coarser lattices, where the effects of the omitted contributions will be larger. 

The construction of the expansion parameter ay is somewhat more complicated, but 
turns out to be equivalent to the following [26|. One first takes the RHS of Eq. (j3]), with 
jj? = e~^/^g^, to define an intermediate coupling, ay(q'^). The corresponding /5 function, 
/?', is then determined to 4-loops by Pq,---,^^. The higher order coefficients, P'^^..., 
however, depend on the presently unknown /?4,5,..., are hence are themselves unknown. 
The final HPQCD coupling, ay, is obtained from a'y by adding terms of 0(af) and 
higher with coefficients chosen to make /9]^ = /?|' = ■ ■ ■ = 0. Since ^4^,5,... are not 
known, the values of the coefficients needed to implement these constraints are also 
not known. The coupling is nonetheless, in principle, well-defined, with higher order 
coefficients computable as soon as the corresponding higher order /^^ become available. 
Since the 4-loop-truncated P^ function is, by defnition, exact, the distortions of the fit 
parameters induced, in general, by the 4-loop truncation of the running are absent for 
the ay coupling. The price to be paid for this advantage is the unknown perturbative 
uncertainty in the relation between ay and as, which affects the conversion and running 
to as{Mz). With this definition, ay differs from ar beginning at 0{a'^). 

The other difference between the two re- analyses lies in the treatment of ri/a, ri, and 
{asC^/n). In Ref. |l|], these are allowed to vary independently, though within the range 
of the input prior constraints, for each Ok, whereas in our analysis, they are treated 
as fixed external input, and have the same central values for all Ok- As noted above, 
the difference in the treatment of ri/a and ri is expected to have a negligible impact. 
The impact of the differing treatments of (a^G^/yr) should be similarly negligible for 
observables with intrinsic scales high enough that the associated correction is small. 

The two different implementations of the original HPQCD/UKQCD approach will 
thus, when restricted to high-scale observables, correspond to isolating residual higher 
order perturbative uncertainties in different sectors of the analysis. If these uncertainties 
are, as desired, small in both cases, the two analyses should be in good agreement. Such 
agreement (which is, in fact, observed, provided comparison is made to the very recent 
HPQCD update) serves to increase confidence in the results of both analyses. 



III. DETAILS OF OUR RE-ANALYSIS 

In our analysis, we have calculated the desired Wilson loops using the publicly available 
a ~ 0.09, 0.12, 0.15 and 0.18 fm MILC nj = 2+1 ensembles and incorporated information 
on Wu and Wu for the three a ~ 0.06 fm USQCD ensembles provided to us by Doug 
Toussaint of the collaboration. 

We follow the basic strategy of the earlier analysis, using the same 3-loop perturba- 
tive input, but with the following differences in implementation. First, we employ the 
expansion parameter ax throughout. All running of a^ is carried out using exact 4- 
loop-truncated running of the intermediate variable, as, whose relation to ay is given by 
Eq. ([2]). Second, to minimize the effect of our incomplete knowledge of the running of ax 
beyond 4-loop order, the impact of which will be larger for coarser lattices, we perform 
"central" 3-fold versions of our fits using the three finest lattices, with a ~ 0.12, 0.09 and 
0.06 fm. Expanded 5-fold fits then serve as a way of studying the impact of the truncated 
running, as well as of the truncation of the perturbative expansion for the Ok- Since we 
do not currently have access to the actual a ~ 0.06 fm configurations, we are restricted 
to analyzing the three observables indicated above. One of these, \og{Wi2/uQ), has a 
significantly lower BLM scale, and hence is particularly useful for studying the impact of 
these truncations. As in Ref. [l|, we extrapolate linearly in the quark masses |27|, and 
estimate (and subtract) residual mass-independent NP effects using the known form of 
the leading order gluon condensate contributions to the relevant Wilson loops. 

Regarding the mass extrapolation, the sets of configurations for different mass combi- 
nations ami /arris corresponding to approximately the same lattice spacing a in fact have 
slightly different measured ri/a. Since the Ok we study are themselves scale-dependent, 
full consistency requires converting the results corresponding to the different am^/ arris to 
a common scale before extrapolation. This could be done with high accuracy if the pa- 
rameters appearing in the perturbative expansion of the Ok were already known. Since, 
however, some of these parameters are to be determined as part of the fit, the extrapola- 
tion and fitting procedure must be iterated. With sensible starting points, convergence 
is achieved in a few iterations. The dominant uncertainty in the converged iterated ex- 
trapolated values is that associated with the uncertainties in ri/a. There is also a 100%- 
correlated global scale uncertainty associated with that on ri. We employ ri = 0.318(7) 
fm, as given in the MILC Lattice 2007 pseudoscalar project update [28|. 

The mass-independent NP subtractions are estimated using the leading order (LO) 



D = 4: gluon condensate contribution, SgWmn, to the m x n Wilson loop, Wmn |19 | 






SgWmn = — — m^n^a^ (a^CVvr) (5) 

3d 

and the central value, {asC^/n) = (0.009 ± 0.007) GeV^, of the updated charmonium 
sum rule analysis [29]. Since the error here is already close to 100%, we take the difference 
between results obtained with and without the related subtraction as a measure of the 
associated uncertainty. This should be sufficiently conservative if the correction is small. 
If not, the measured Ok values may contain additional non- negligible mass-independent 



contributions, of dimension D > 4, which we do not know how to estimate and subtract. 
Ok for which this occurs will thus provide a less reliable determination of ag- 

Fortunately, for the observables we consider, the gluon condensate correction is, as 
desired, small. For Ok = log(Wii), the corrections required for the 3-fold (5-fold) fit do 
not exceed ~ 0.1% (~ 0.5%). The corrections remain small (less than ~ 0.4% (~ 1.8%)) 
for Ok = log(iyi2). The effect is somewhat larger for \og{Wi2/uQ), as a consequence of 
cancellations encountered in combining the uncorrected logfiyn) and log(iyi2) values, 
but still reaches only ~ 1.3% (~ 5.6%) for the 3-fold (5-fold) fit H- 

In line with what was seen in Ref. [l|, we find that the known terms in the per- 
turbative expansions of the Ok are insufficient to provide a description of the observed 

(k) 

scale-dependence, even when only the three finest lattices are considered. When Cg is 
added to the fit, however, we find very good fits, with x^/dof < 1 (very significantly so 
for the 3-fold fits). With current errors, it is thus not possible to sensibly fit additional 
coefficients c^>3. This raises concerns about possible truncation uncertainties. Compar- 
ison of the results of the 3-fold and 5-fold fits provides one handle on such an uncertainty 
since the relative weight of higher order to lower order terms grows with decreasing scale. 
If neglected higher order terms are in fact not negligible, the growth with decreasing scale 
of the resulting fractional error should show up as an instability in the values of the pa- 
rameters extracted using the different fits. We see no signs for such an instability within 
the errors of our fits, but nonetheless include the difference of central values obtained 
from the 3-fold and 5-fold fits as a component of our error estimate. 



IV. RESULTS 

Central inputs for our fits are the measured lattice observables (whose errors are tiny 
on the scale of the other uncertainties), the computed Dk, q and ^ [l|,ll6|, ri/a, ri and 
{asC^/n), and the choice of the 3-fold fitting procedure. In addition to the uncertainties 
generated by the errors on ri/a, ri and (agG^/vr), are those due to uncertainties in 
numerical evaluations of the Dk , c^ and Cg . 

We construct an "overall scale uncertainty error" by adding linearly the fit uncertain- 
ties generated by those on ri and the ri/a. This combined error is added in quadrature 
to (1) uncertainties produced by varying the C2 (and, if relevant, c{ ) within their er- 
rors, (2) the difference between results obtained with and without the gluon condensate 
correction, and (3) the difference between the results of the 3-fold and 5-fold fits. Be- 
cause of the iterative nature of the fit procedure, the mass extrapolation uncertainty is 
incorporated into what we have here identified as the overall scale uncertainty. 

We run our nf = 3 results to Mz using the self-consistent combination of 4-loop 
running and 3-loop matching at the fiavor thresholds, taking the fiavor thresholds to 
lie at rmdrnc) and rmb{rnb), with mc{mc) = 1.286(13) GeV and mb{mb) = 4.164(25) 



GeV [31|, and r allowed to vary between 1 and 3. These uncertainties in the matching 
thresholds, together with standard estimates for the impact of the truncated running and 
matching, produce an evolution contribution to the uncertainty on as{Mz) of ±0.0003 



TABLE II: Central fit results for as{Mz) and the Cg 



Ok 


as{Mz) 


^3 


log(PFii) 

log (Wu) 

log {W12/4) 


0.1192 ±0.0011 
0.1193 ±0.0011 
0.1193 ±0.0011 


-3.8 ±0.6 
-4.0 ±0.9 
-1.7 ±0.8 



Our central fit results for as{Mz) and the c^ are given in Table HTl For comparison, 
the results for a,(Mz) obtained in Ref. pL] were 0.1171(12), 0.1170(12) and 0.1162(12), for 
\og{Wu), logiWu) and log{Wi2/uQ), respectively. Our a,^(Mz) are significantly larger, 
and in closer mutual agreement. The recent HPQCD update 17| also finds significantly 
larger values. (We will return to a more detailed comparison of the two updates below.) 
The very good agreement between the as{Mz) values obtained in our fits using both low- 
and high-scale observables suggests that the effects of the truncated running, present at 
some level in all such fits, are small in the cases we have studied. 

One-sided versions of the various components of the total errors on as{Mz) are dis- 
played in Figure [H The difference of the 3-fold and 5-fold determinations is ~ 0.0004, 
significantly smaller than the ~ 0.0009 overall scale uncertainty. The results thus show 
no evidence for any instability associated with opening up the fit to lower scales. 

While the total error on as{Mz) is the same for all three Ok considered, the general 
arguments above lead us to believe that the most reliable determination is that obtained 
using the highest-scale observable, log(W^ii), and highest-scale (3-fold fit) analysis win- 
dow. Our final assessment, 



as{Mz 



0.1192 ±0.0011 



(6) 



is in excellent agreement with the non-lattice average and the result, 0.1184 ± 0.0009, 
of the independent HPQCD analysis. The various results are shown for comparison in 
Figure [21 A more detailed discussion of the relation between our re-analysis and that of 
HPQCD may be found in the Appendix. 
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FIG. 1: Contributions to the errors on 
as{Mz)- Shown are the as{Mz) obtained us- 
ing (i) the 3-fold fit strategy, with all cen- 
tral input, (ii) the alternate 5-fold fit strat- 
egy, with all central input, and (iii) the 3-fold 
fit strategy, with, one at a time, each input 
shifted from its central value by la, retaining 
central values for the remaining input param- 
eters. The error bars shown are those associ- 
ated with the uncertainties in ri/a. 
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cent non-lattice determinations. 



APPENDIX: MORE ON THE RELATION TO THE HPQCD RE- ANALYSIS 

After the completion of the work reported in this paper, the HPQCD Collaboration 
posted an update of their earlier 2005 analysis [l7|. This update works with a subset of 
11 of the available MILC ensembles, spanning the a ~ 0.18, 0.15, 0.12, 0.09 and 0.06 fm 
lattices and a range of ami/ arris. The fits follow the strategy of the earlier analysis pLJ], 
employing the expansion parameter ay, and fitting the unknown Cn using priors. Linear 
mass extrapolation has been employed, and mass-independent NP D = A contributions 
estimated and subtracted using the LO formula for SgWmn- The fitting of ri/a, ri and 
{agC^/ii), observable by observable, using central input and prior widths, is also as in 
the earlier analysis, with the exception that the central value and width for (asG^/vr) 
are now and ±0.012 GeV^, respectively. 

The HPQCD implementation differs from ours in the choice of expansion parameter, 
and in the implementation of the input information on ri/a, ri and (— G^). For the 
reasons discussed above, we expect the impact on as{Mz) of the observable- by-observable 
fitting of ri/a, ri and (asG^/ir) in the HPQCD approach to be small for Ok having 
small gluon condensate corrections. Since the different choices of expansion parameter 
correspond to different ways of isolating residual higher-order perturbative uncertainties. 
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one expects the results of the two analyses to be in good agreement so long as (i) one 
is working with Ok having small mass-independent NP corrections, (ii) the same input 
values are used for both, and (iii) residual NP and higher-order perturbative uncertainties 
are indeed small. The situation is likely to be more complicated for Ok with sizeable 
estimated D = A gluon condensate corrections. 

The results of the HPQCD fit for the three Ok we consider are as{Mz) = 0.1186(9), 
0.1186(9) and 0.1183(8) for log(iyii), log(iyi2) and \og{Wi2/ul), respectively 0. AH 
are in good agreement within errors with the corresponding results from our analysis. 
This agreement is further improved if one takes into account the small difference in 
input Ti values. Were we to switch from ri = 0.318 fm to the central value of the 
HPQCD determination, 0.321(5) fm, all three of our as{Mz) results would decrease by 
0.0002. Note also that use of the central charmonium sum rule input for (asC^/ir) in our 
calculation raises the output as{Mz) obtained from \og(Wu), logiW^) and logiWu/ul) 
by 0.0001, 0.0004 and 0.0005, respectively. Our fitted values would thus be in even closer 
agreement with those of the HPQCD update were we to impose the HPQCD central zero 
value of (ttsG^/vr) in our fits. Since the fitted (agC^/n) values obtained by analyzing the 
various Ok are not quoted in Ref. [l7J, it is not possible to quantify further the role of 
this effect. The agreement for the three observables under discussion is in any case good, 
within expectations, independent of this question. 

We now turn to a more detailed discussion of the issue of the subtraction of the 
mass- independent NP contributions. If the estimated LO, D = 4 gluon condensate 
subtraction represents only a small fraction of the measured Ok at the scales under 
consideration, analogous mass-independent NP contributions with D > A should be even 
smaller, and hence safely negligible. If, however, the estimated D = A correction is 
sizeable, analogous D > 4 corrections can no longer be expected to be small. These 
necessarily scale differently with lattice spacing than do the D = perturbative and 
D = 4 NP contributions and hence, if not included when fitting the data, are likely to 
force shifts in both a^ and {agG^/n) if present at a non-negligible level. 

We deal with this potential problem by focussing on Ok for which the impact of the 
estimated D = 4 gluon condensate subtraction is small compared to the variation of 
the Ok in question over the lattice scales employed in the fit. In the initial version of 
the HPQCD re-analysis, mass-independent NP subtractions were estimated using only 
the D = 4 gluon condensate form, even for observables where the estimated correction is 
sizeable. In the more recent update, additional terms, scaling as would mass-independent 
contributions oi D > 4, are added to the fit function for each observable, and the ac- 
companying coefficients extracted as part of the augmented Bayesian fit. The impact of 
including the D > 4 terms is, as expected, small for those observables having small values 
of the estimated D = 4 subtraction. For observables with larger D = 4 subtractions, the 
fit errors are increased (by factors of ~ 2 for those observables having the largest D = 4 
corrections) and some shifts in as{Mz) of order 1/2 to 1 times the smaller preliminary 
errors are observed. The shifts serve to reduce the spread of as{Mz) values compared 
to that seen in the original version of the re-analysis. The values of (asC^/n) obtained 
from the independent fits to the different observables are not quoted in Ref. [17|], but a 
useful test of the self-consistency of the approach would be to verify that the inclusion of 
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the D > 4 contributions has brought these values into good agreement with one another. 

It is worth noting that the observables, \og(W23/ ul^), log(Wi4/W23), and 
log(iyii 1^23/^^12^^13)5 which produce the three smallest results for as{Mz), have esti- 
mated D = 4 corrections significantly larger than those for any of the other observables. 
The magnitudes of the corrections in these cases represent ~ 50 — 100% of the variation 
with scale of the uncorrected Ok between the lightest mass a ~ 0.06 and a ~ 0.12 fm 
ensembles. (This variation-with-scale provides a suitable measure for use in assessing the 
importance of NP corrections since it is the variation with scale which provides the input 
needed to fix the fit parameters, and, as explained in Ref. [l7|, the a ~ 0.06, 0.09 and 0.12 
fm ensembles which dominate the HPQCD re-analysis.) These observables are thus, for 
the purposes of the analysis, rather non-perturbative. Were one to exclude observables 
with larger NP contributions from the HPQCD average, on the grounds that the related 
subtractions introduce additional theoretical systematic uncertainties, the HPQCD re- 
sult would be brought into even closer agreement with ours, though the resulting shift 
would in fact be small (at the ~ quarter a level). 

We stress that, independent of these questions, our results agree well within errors with 
those of the HPQCD update. This agreement is further improved by a shift to common 
input. We argue that the non-zero central value for (a^G^/Tr) obtained from the updated 
charmonium sum rule analysis represents our best present knowledge of this quantity, 
and hence also the best choice as input for evaluating the small mass-independent NP 
subtractions needed for extracting as{Mz). In addition, for the reasons just discussed, 
we believe that the most reliable determinations of as{Mz) are those based on those 
observables for which the D = 4 correction is as small as possible. Such an assessment 
produces the results already noted above, which are in extremely good agreement with 
what it known from other sources. 
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